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Q '. ABSTRACT 

in ■ 

(y>j Aims. We present a wide-field multi-color survey of a homogeneous sample of eleven clusters of galaxies for which we measure total masses 

and mass distributions from weak lensing. 
^ Methods. The eleven clusters in our sample are all X-ray luminous and span a narrow redshift range at z = 0.21 ± 0.04. The weak lensing 
• i-H , analysis of the sample is based on ground-based wide-field imaging obtained with the CFH12k camera on CFHT We use the methodology 

■ developed and applied previously on the massive cluster Abell 1689. A Bayesian method, implemented in the Im2shape software, is used to fit 
5— ( the shape parameters of the faint background galaxies and to correct for PSF smearing. A multi-color selection of the background galaxies 

is applied to retrieve the weak lensing signal, resulting in a background density of sources of ~ 10 galaxies per square arc minute. With the 
present data, shear profiles are measured in all clusters out to at least 2 Mpc (more than 15' from the center) with high confidence. The radial 
shear profiles are fitted with diiferent parametric mass profiles and the virial mass M200 is estimated for each cluster and then compared to 
other physical properties. 

Results. Scaling relations between mass and optical luminosity indicate an increase of the M/L ratio with luminosity (M/L oc L"**) and 
a Lx - M200 relation scaling as Lx 0^ ^200'*" " while the normalization of the M200 relation is close to the one expected from 

hydrodynamical simulations of cluster formation as well as previous X-ray analyses. We suggest that the dispersion in the M200 - and 
M200 - Lx relations reflects the different merging and dynamical histories for clusters of similar X-ray luminosities and intrinsic variations in 
their measured masses. Improved statistics of clusters over a wider mass range are required for a better control of the intrinsic scatter in scaling 
relations. 
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1. Introduction 

Clusters of galaxies are potentially powerful probes for cos- 
mology. They form the high-mass end of the mass function of 
collapsed halos, whose development as a function of redshift is 
a basic test of the hierarchical structure formation scenario and 
depends sensitively on a nu mber of cosmological parameters 
dEke et al. I ll996t IVoitll2005h The distribution of mass within 
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clusters fo rms another test of the non-hnear development of 
structures (INavarro et al.lll997h . The main difficulty in the ap- 
plication of these tests is to accurately measure total masses and 
the distribution of mass in clusters from observational proxies 
which are only more or less indirectly related to mass or the 
gravitational potential. 

Several observational techniques are available to probe the 
mass distribution in clusters, each of them based on different 
physical assumptions and having its own strengths and weak- 
nesses. Observation of the internal dynamics of clusters, based 
on the virial theorem and using the cluster galaxies as test parti- 
cles of the cluster potential, is the "historical" approach, which 
provided early evidence for th e existence of "missing" (now 
"dark") matter dZwickv 19371) . However, clusters of galaxies 
are far from being simple relaxed systems and their structural 
complexity makes analysis of the vel ocity field difficult , as 
soon as the system shows substructure (ICzoske et al ]l2002h . 

Due to its dependence on the square of the electron den- 
sity. X-ray emission from the hot intra-cluster gas (IGM) traces 
the deeper parts of the cluster potential and can be used to in- 
fer the total cluster mass under the assumption of hydrostatic 
equilibrium. It is known, however, that hydrostatic equilib- 
rium alone provides an incomplete description of the physics of 
the IGM. X-ray observations with XMM-Newton and Chandra 
have rev ealed a wealth of complexity in the cluster X-ray 
emission ( Markevitch & Vikhlininl2007i) which sign-post devi- 
ations from hydrodynamic equilibrium on various levels rang- 
ing from "sloshing" in apparently r elaxed clusters to gross 



2007t ICorless & Kingl l2007h to large errors in the mass es- 
timates in cases where there are unrec ognized line-of-sigh t 



deviations in, e.g ., cluster mergers (Finoguenov et al.l 12005 



Clowe et al.ll2006l) . iNagai et al.l (120071) show from simulations 



that mass estimates of relaxed clusters are biased low by 5 to 
10% under the assumption that the gas is supported only by hy- 
drostatic pressure; for unrelaxed clusters the situation is much 
worse. Hence, the usual X-ray observables like Lx or Tx do 
not provide the simple and expected robust mass estimators, 
although the new ly introduced Fx para meter (i.e. the product 
of Tx and M„ snn. lKravtsov et al.l (120061) ) shows some promise, 
with a low scatter when scaled with the mass. 

The measurement of the Sunyaev-Zeldovich (SZ) eff'ect 
is sensitive to the integrated pressure of the intra-cluster 



gas and is potentially a ro bust mass estimator (Greg o et al 



2001 : Bonamente et al ■|2006 ). The new SZ interferometers like 
AMI dBarker et all l2006l AMI collaboration) and millimeter 
bolometer arrays (e.g. LABOCA on APEX) are promising fa- 
cilities for measuring the SZ effect in clusters, although their 
limited spatial resolution will limit the quality of cluster mass 
distribution. 

Finally, weak gravitational lensing is most directly related 
to the gravitational potential and hence the total mass distri- 
bution; in particular, lensing does not rely on any assumptions 
concerning the physical state of the system. The main system- 
atic effect affiicting lens mass measurements is that lensing 
measures projected masses, hence the interpretation of lensing 
measurements in terms of physical three-dimensional masses 
relies on certain assumptions concerning the spatial distribution 
of the matter Projection effects can range from slight biases 
arising from triaxiality of the clusters or projection of the gen- 
eral large-scale structure (IMetzler et aLluOOlt iKing & Corless 



mergers of clusters of comparable size (ICzoske et al.ll2002h . 



Furthermore, there are calibration issues in the mass measure- 
ments coming from 1) the removal of instrumental effects act- 
ing on the galaxy shapes, 2) the calibration of the weak lensing 
signal, 3) the dilution of the signal by contamination with faint 
cluster or foreground galaxies, and 4) the uncertainties in the 
source redshift distribution. The first two issues have largely 
been addressed with simulated datajnjhe STEP collaboration 
(Shear TEsting Program, iHevmans et al.ll2006l) aimed at com- 
paring many different weak lensing methodologies and iden- 
tifying systematic effects in the respective methods. The last 
two issues can be addressed if additional information is avail- 
able such as multi-color photometry and photometric redshifts 
or spectroscopic data. 

From the above discussion it is clear that joint analyses of 
several types of observations are necessary if one wants to fully 
understand the structure of clusters of galaxies and the relation 



between different observab les and mass estimates (IDahle et al 



20021: ICypriano et al ]|2004 . In practice, such an in-depth anal 



ysis can only be conducted for comparatively small samples of 
clusters, from which the relations and the scatter around the 
relations have to be empirically calibrated. 

The aim of our program is to study a homogeneous sam- 
ple of galaxy clusters using a variety of te chniques to constrain 
their mass distribution. ISmithet al.l(l2005b presented the results 
from strong lensing mass modeling of HST/WFPC2 observa- 
tions of the central parts of the clusters and a comparison to 
X-ray temperatures measured by Chandra. Analyses of XMM- 
Newton observations of mostly the same clusters are presented 
by Zhang et al. (2007). In this paper, we present the weak lens- 
ing analysis of the cluster sample, usin g imaging data obtained 



with the CFH12k wide-field camera dCuillandre et al.1 l2000l) 



mounted at the Canada-France-Hawaii Telescope (CFHT); 
eleven clusters were observed in the B, R and / bands. The 
weak lensing analysis of the full sample of clusters follows 
the methodology presented previously by Bardeau et al. ( 20051 
hereafter Paper I) who also applied it to the cluster A 1689. 

This paper is organized as follows: Sect. |2] presents the 
cluster sample and gives a summary of the data reduction and 
the construction of catalogs that are used in the weak lensing 
and optical analyses of the clusters. In Sect. [3] we convert the 
galaxy shape measurements into two-dimensional mass maps 
which allow a qualitative assessment of the morphology of the 
dark matter distribution and comparison to the galaxy distri- 
bution in the clusters. In Sect. |4] we provide quantitative mea- 
surements of the cluster masses using one-dimensional radial 
fits to the weak shear signal. Sect. |5] discusses the light distri- 
bution of the clusters compared to the weak lensing mass in 
order to measure the mass-to-light ratio. We also compare the 
weak lensing masses to the X-ray properties of the clusters (lu- 
minosity and temperature) and normalize the M-Lx and M-Tx 
relations. We summarize our conclusions in Sect.|6]and discuss 
prospects for future weak lensing cluster surveys. Finally we 
give some details on the individual properties of the clusters in 
the appendix. 
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Throughout the paper we use Ho - 70 km s ' Mpc ' , = 
0.3, = 0.7. At a redshift of z = 0.2, 1" corresponds to 
3.3 kpc and 1' to 200 kpc. Magnitudes are given in the Vega 
system. 

2. Observations and data reduction 

2.1. Description of the cluster sample 

The clusters anal yzed in this paper w ere selected from the 



XBACs catalog of lEbeling et alj (1199 6). a flux-limited compi 



lation of Abell clusters detected in the Rosat All-Sky Survey 
data. While th e sample is thus b ased on the optically selected 
Abell catalog ( Abell etaPl 19891) . which is known to be incom- 
plete at high r edshift and lo\y mass, comp arison with the X-ray 
selected BCS dEbeling et al.lll998il2000l) shows that more than 
74% of all BCS clusters are indeed Abell clusters. Since, for 
the very X-ray luminous systems considered here, this fraction 
rises to almost 90%, the optical preselection should not bias 
our sample. We select systems within the relatively narrow red- 
shift range of 0.17 < z < 0.26 in order to obtain an approxi- 
mately luminosity-limited sample and to ensure that all clus- 
ters lie at about the same distance from the background popu- 
lation. Si nce X-ray luminosity is broa dly correlated with clus- 
ter mass (IReiprich & Bohringerll2002h . our sample is approxi- 
mately mass limited, too. The redshift range quoted above was 
chosen to maximize the lensi ng efficiency for a backg round 
galaxy population at (z) ~ 1.0 ( Nataraian & Kneiblll997 ). thus 
defining the requirements in terms of limiting magnitudes. We 
applied further limits in declination (-20° < 6 < 60° for 
accessibility with CFHT from Mauna Kea), Galactic latitude 
(\b\ > 20° to minimize contamination by stars) and hydrogen 
column density (A^h < 10 x 10^° cm"^) and from the remaining 
clusters randomly selected a sample of twelve clusters within a 
redshift dispersion cr./z = 12%. Of these, eleven clusters were 
observed with the CFH12k wide-field camera (A773 could not 
be observed within the allocated nights). The same cluster sam- 
ple was also imaged with the Hubble Space Telescope (HST) 
in order to measure the mass distribution in the cluster cores 
using strong-lensing techniques. Eight systems were observed 
by us with the WFPC2 camera (Program ID: 8249, PI Kneib), 
and data from obser vations of another two were retrieved from 
the HST archive (see lSmith et al .l2005h . Note that sample com- 
pleteness is not of critical importance for this project which 
rather aims at compiling a homogeneous data set for a represen- 
tative comparison among mass measurement techniques. Fig.II] 
shows how our sample covers the high-luminosity region in the 
XBACs catalog. 

Table [T] provides a summary of the cluster properties in- 
cluding global X-ray characteristics that will be used in our 
analysis. 



2.2. Observations 

Imaging data were obtained at the Canada-France-Hawaii 
Telescope during three observing runs in February 1999, 
November 1999 and May /June 2000 using the CFH12k cam- 
era, a mosaic of twelve 2k x 4k CCDs dCuillandre et al.l2000,) . 




A Sample member 
X > 10 

+ 6 > 60° or d < -20° 

* \b\ < 30° 



Fig. 1. Lx-z diagram for clusters from the XBACs catalog 
(dots). Triangles mark the 12 members of our sample, although 
one cluster (A 773) could not be observed and is not studied 
in this paper The other symbols mark clusters which failed 
some of our secondary selection criteria. The solid line corre- 
sponds to an X-ray flux of 5 x 10"'^ erg s"' cm"^, the flux limit 
of XBACs. The point outside the redshift limits (dashed lines) 
corresponds to the cluster A2218. Luminosities given here are 
the original XBACs luminosities from Rosat All-Sky Survey 
data computed for Ho = 50kms"' Mpc"', Qm = 1, i^A - 0. 



With a pixel scale of 0'.'206 this camera covers a field of view of 
42 X 28 arcmin^, i.e. about 1/3 of a square degree. Observations 
were conducted in three filters, B, R and /. A1835 was observed 
in V instead of B in February 1999 because the latter filter was 
not yet available at the time. Further observational details are 
given in Table |2] 

We determine the limiting magnitude of each final image by 
computing the magnitude of a point source detected at 5cr in an 
aperture with diameter 1 .45 times the full width at half maxi- 
mum (FWHM) of the seeing disk. As shown by the CFH12k 
exposure time calculator (DIElQ) this is the optimal aperture 
within which 96% of the flux is integrated while the noise level 
remains sufficiently low. A point-like object detected at this 
significance level will have a magnitude error of about 20%. 
The resulting limiting magnitudes are included in Table |2] The 
measured magnitude limits corresponding to 50% complete- 
ness are found to be about 1.8 magnitudes brighter, partly due 
to the small apertures used in DIET. 

During all three runs, the observing conditions were very 
good in terms of seeing (a prime requirement for weak-lensing 
studies) with the FWHM of stars consistently below 1" in the 
R band. The first and second observing runs were photometric; 
the third run was affected by cirrus (see Sect. l2.3l i. 

2.3. Data reduction 



Full d etails of the data reduction are presented in ICzoske 
(l2002l) . n Paper I we give a summary of the reduction process 
with special emphasis on the astrometric alignment of the ex- 
posures, which is of particular importance for the weak-lensing 



|http : //www ■ cfht . hawaii . edu/Instrxments/Imaging/CFH12K/DIET/CFH 
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Table 1. Physical properties of the eleven clusters of the sample studied in this paper The coordinates are those o f the central 
galaxy as measured on the CFH12k images. The X-ray luminosities are taken from XMM-Newton measurements in lZhang et al. 
(120071) and are X-ray bolometric luminosities excluding the < O.lrsoo region. The X-ray temperatures listed in column 6 were 
measured from Chandra data by Smith et al. (2005), excluding the very central regions where cool cores might be present in 
some of the clusters, except for A1689* dAndersson & Madeisldil2004 XMM data^ and A2390** dAllen et al.ll200li) . Column 7 
lists the X-ray temperature measured with XMM-Newton from Zhang et al.l (l2007h . except for A 2219'** for which we list the 
temperature obtained with ASCA (Ota & Mitsuda 2004) as no XMM data are available. 



Cluster 


RA 
(J2000) 


Dec 
(J2000) 




z 


Lx (XMM) 

W^h-l ergs-' 


Tx (Chandra) 
keV 


Tx (XMM) 
keV 


A 68 


00'"37™0ff.9 


-^09°09'24" 





255 


10.1 ±0.9 


9-5 !|:^ 


7.7 ±0.3 


A 209 


01''31"52'!6 


-13°36'40" 





206 


I3.2± 1.1 


8.7 


7.1 ±0.3 


A 267 


0li>52™42W 


+01°00'26" 





230 


6.6 ±0.7 


6 


6.5 ±0.4 


A 383 


02''48"03!4 


-03°31'45" 





187 


4.6 ±0.5 


5 2 

-'■^ -0.2 


5.3 ±0.2 


A 963 


10''17"03!6 


+39°02'50" 





206 


10.2 ±0.9 


7 9 +0.3 
' -0.3 


6.3 ± 0.2 


A 1689 


13''11"30!1 


-01°20'28" 





184 


21.4± 1.0 


9 +013 • 
-0.12 


8.4 ±0.2 


A 1763 


13''35"20!1 


+4r00'04" 





228 


15.9 ± 1.4 


7 7 +0.4 
'■' -0.4 


6.3 ±0.3 


A 1835 


14'"01"'02':l 


-t-02°52'42" 





253 


30.0 ± 1.4 


Q ri +0.6 
y-^ -0.4 


8.0 ±0.3 


A22I8 


16''35™5H5 


+66° 12' 15" 





171 


12.2 ±0.9 


6-8 


7.4 ±0.3 


A2219 


16M0"19i9 


+46°42'41" 





228 




13.8 


9.2 ± 0.4 


A 2390 


21''53"'3ffi9 


+ 17°41'43" 





233 


28.9 ±2.2 


11.5 !li " 


10.6 ± 0.6 



Table 2. Summary of the CFHT observations of our cluster sample. Exposure times are given in seconds. The seeing is measured 
as the FWHM of stars in the final stacked image for each cluster. The limiting magnitudes were computed following the DIET 
recipes (see text for details) for a point source detected at 5cr in an area of diameter 1.45 times the FWHM of the image. At this 
limit magnitudes are measured with an intrinsic error Am ~ 0.2. They do not include the correction for galactic absorption which 
is given for each cluster in the 3 filters. 



Cluster 


Integr. 
time 


Seeing 


B 

Limiting 
mag. 


Ab 


Integr. 
time 


R 

Seeing 


Limiting 
mag. 




Integr. 
time 


I 

Seeing 


Limiting 
mag. 


A, 


A 68 


8100 


I'.'l 


26.4 


0.40 


7200 


a' 7 


26.3 


0.25 


3600 


a' 6 


25.3 


0.18 


A 209 


7200 


I'.'O 


26.8 


0.08 


6600 


a' 7 


26.3 


0.05 


3600 


a' 7 


25.1 


0.04 


A 267 


3000 


I'.'O 


26.0 


O.ll 


4800 


a:i 


25.9 


0.07 


900 


a' 7 


24.5 


0.05 


A383 


7200 


a' 9 


27.0 


0.14 


6000 


a' 9 


26A 


0.09 


3600 


a' 7 


25.1 


0.06 


A 963 


7200 


a' 9 


27.0 


0.06 


4800 


a' 8 


26.1 


0.04 


10500 


I'.'l 


24.6 


0.03 


A 1689 


3600 


a' 9 


26.7 


0.12 


3000 


a' 8 


26 A 


0.07 


3000 


a' 9 


24.8 


0.05 


A 1763 


3600 


I'.'O 


26.7 


0.04 


6000 


a' 9 


26.2 


0.02 


3000 


a' 8 


25.0 


0.02 


A 1835 


3750* 


a.' 8 


26.4* 


0.10* 


5400 


a' 7 


26.5 


0.08 


3750 


a' 8 


25.5 


0.06 


A22I8 


3378 


I'.'l 


26.3 


0.11 


6900 


I'.'O 


26.2 


0.07 


3000 


a' 8 


24.7 


0.05 


A22I9 


5400 


I'.'O 


26.8 


0.11 


6300 


^.'8 


26.4 


0.07 


3000 


a' 8 


25.1 


0.05 


A 2390 


2700 


I'.'l 


26.2 


0.48 


5700 


a' 7 


26.3 


0.30 


3600 


a: 9 


25.1 


0.22 



* V-band data for A 1835. 



analysis of the images. We here provide additional information 
on the absolute photometric calibration of the images which 
will be relevant for the comparison of the total galaxy cluster 
luminosity to the total mass distribution presented in Sect. |5] 

Before combining the exposures taken through the same fil- 
ter for a given cluster we apply empirically determined photo- 
metric scaling factors to the exposures, which take into account 
differences in air mass and atmospheric transparency. The fac- 
tors are determined by comparing the instrumental fluxes of 
several thousand stars in each image to the corresponding 
fluxes in a reference exposure (usually the exposure taken at the 



lowest air mass). The median of the distribution of flux ratios is 
selected and used as the global scale factor which brings the ex- 
posure to the same level as the reference exposure. Remaining 
variations in the photometric zero point between chips (which 
have not been completely corrected by the twilight flat fields) 
are determined in a simpler manner using the sky levels in ad- 
jacent parts of neighboring chips. The intrinsic accuracy of the 
photometric calibration within the images is estimated to be on 
the order of 0.01 mag. 

Plotting the relative scaling factors for the individual expo- 
sures against air mass allows the determination of the atmo- 
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spheric extinction if the series of exposures has been taken 
over a sufficient range in air mass; this is the case for sev- 
eral of our fields. We confirm that for the first and sec- 
ond of our observing runs the conditions were photomet- 
ric throughout; the extinction coefficients determined in this 
way are in good agreement with the values listed on the 
CFH12k web siteB The photometric zero points are deter- 
mined from observations of photometric stand ard field s taken 
throughout the nights. We use th e catalogs by Landolti (119921) 
supplemented more recently by IStetson ( 200C ). Finally, the 
Ga lactic extinctio n is co mputed for each cluster from the maps 
of ISchlegel et aP (fl998). accessed through the NASA/IPAC 
Extragalactic Database Extinction corrections in R 

are typically on the order of 0.05 mag with the exception of 
Abell68 and Abell2390 for which we find values of up to 
0.30 mag. 

The third observing run was affected by intermittent thin 
cirrus and the conditions were not photometric. Calibrating 
the images taken during this run with standard observations 
is therefore likely to result in systematic errors in the photo- 
metric zero points. For the purposes of the present paper we 
adjust the zero points a posteriori by selecting stars from our 
fields and comparing their col or distribution to external multi- 
color photometric sequences jSaha et al. I l2005b . In the 7? - / 
\s. B - R color plot the stellar sequence displays a character- 
istic knee from which photometric color corrections in B - R 
and R - I can be determined. A correction of 0.1 to 0.4 magni- 
tudes is required for the clusters A1689, A1763, A2218, A2219 
and A2390 to obtain a correct value for the R - I color of el- 
liptical galaxies. The correction terms are arbitrarily applied 
relative to the R-band which is taken as reference. Note that 
except for the color selection of cluster ellipticals or of "red" 
galaxies (see below for their definition), only the R-band lumi- 
nosities, which are only mildly affected by the non-photometric 
conditions, are used in the rest of this paper Only the clusters 
A 1763 and A 2219 suffered from significant absorption during 
the R-band observations, with an unknown absorption factor 
not higher than 0.3 mag. 

2.4. Catalogs 

From the corrected images, photometric catalogs are created 
for each cluster using SExtractor ( Bertin & Arnouts 1996). 
Different catalogs are produced for the various measurements 
needed for the analysis. The details of the general methodology 
are given in Paper I and are summarized below. Object colors 
are computed using aperture magnitudes in circles of 3" diam- 
eter, and stellar objects are identified using the tight relation in 
the magnitude vs. peak surface-brightness diagram. 

Color-magnitude diagrams are constructed from the photo- 
metric catalog (Fig. |2]i, and the cluster red sequence is manu- 
ally identified and approximated by a line with negative slope; 
all galaxies in a certain range around the line (Am = ±0.08 
typically) are stored in a catalog which we will refer to as the 
"elliptical galaxy" catalog. Since, for simplicity's sake, we do 



A 1763 




Fig. 2. Color-magnitude diagram for the cluster Abell 1763 
showing the characteristic red sequence of elliptical cluster 
galaxies (width Am - +0.06). All galaxies redder than the se- 
quence are included in the "red" galaxy catalog. The large dots 
correspond to galaxies within a radius of 250" while the small 
dots come from the full catalog. 

Table 3. Photometric properties of the sub-catalogs used in the 
weak-lensing analysis. Only R-band data are summarized. The 
magnitude range used for the faint galaxy selection is defined 
to be: [m^ + 3.5 to mco^p + 0.5], nJcomp being the completeness 
magnitude of the catalog. 



Cluster 


m* 


Faint galaxies 
magnitude range 


Mean number density 
(gal / arcmin^) 


A 68 


19.04 


[22.5-24.9] 


22 


A 209 


18.41 


[21.9-24.9] 


24 


A 267 


18.71 


[22.2-24.5] 


17 


A 383 


18.12 


[21.6-24.9] 


22 


A 963 


18.41 


[21.9-24.8] 


23 


A 1689 


18.08 


[21.6-24.8] 


22 


A 1763 


18.67 


[22.2-25.0] 


24 


A 1835 


19.02 


[22.5-25.4] 


25 


A2218 


17.91 


[21.4-25.0] 


21 


A2219 


18.67 


[22.2-25.3] 


25 


A 2390 


18.75 


[22.3-24.6] 


18 



not introduce any cut in magnitude, this catalog is affected by 
field contamination, especially at faint magnitudes. 

For each cluster the photometric catalog is split into two 
sub-catalogs containing "bright" and "faint" galaxies, respec- 
tively. The magnitude cuts used to select these catalogs are 
scaled to the value of ;«*, determined separately for each clus- 
ter and each filter (see Paper I for details and Table [3]). The 
"faint galaxy" catalog is dominated by background sources 
and is used for the weak-lensing analysis. The galaxy density 
ranges from 22 galaxies/arcmin^ in 7? to 14 galaxies/ arcmin^ in 
B and /. These numbers are averaged densities over the cluster 
sample and differ slightly from cluster to cluster due to cos- 
2 |http: //w ww, cfht .hawaii . edu/I nstriments/Imaging/CFH12H >lgumaay9^eMrgBFC ^£t^te tMn the photometric depth of the 
^ |http : //n'edwww . ipac . caltech . eduT] catalogs (TableO. 
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Fig. 3. Observed color index R - I versus redshift for differ- 
ent spectral types, ranging from elliptical (with and without 
evolution, marked by full and dotted lines, respectively) to 
the bluest star-forming galaxies. We select "red galaxies" with 
R - I > 0.7 (outside the hatched region): The resulting catalog 
is next to free from foreground galaxies and cluster members, 
and is limited to z 5. 1-8 . Synthetic spectral models are from 
Bruzual & CharlotI [20031) . 



To further select faint background galaxies as well as to 
minimize cluster contamination, we add a simple color cri- 
terion to isolate galaxies above the cluster red sequence. As 
shown in Fig. [3] galaxies redder than R- 1 ^ 0.7 should essen- 
tially be background galaxies ranging from the cluster redshift 
up to z ^ 1.5 - 1.8. This approach has alre ady been success- 
fully developed and use d by several groups jKneib et al ] l2003l 



Broadhurst et al l2005b ^ to eliminate or at least severely reduce 
the contamination of background galaxy samples by cluster 
galaxies. This extra criterion reduces the galaxy number densi- 
ties to typically 8 - lOgalaxies/arcmin^ (a reduction by a fac- 
tor two compared to the "faint galaxy" catalogs, see Fig. |6]l. 
However, because the "red galaxy" catalog can be assumed to 
contain mostly background sources, we use this catalog for the 
remainder of our weak-lensing analysis, except in Sect.[3]where 
for practical reasons the "faint galaxy" catalogs are used as in- 
put. 

2.5. Measurement of galaxy shape parameters 

The shape parameters of faint galaxies are measured in all im- 
ages, allowing a qualitative comparison between the different 
bands. However, in the end we only retain the measurements 
in the R band because the data quality in this band in terms 
of seeing and source density is superior to that of the other 
bands. In order to measure the galaxy shapes for the weak- 
lensing analysis, we correct them fo r PSF sm earing using the 
Im2shape software developed by Br idle et al. (2001). Details of 
the implementation are given bv Hevmans et al.l (12006) as well 
as a comparison wit h more commonly used techniques such 



V 

as the KSB scheme dKaiser et al.lll995l) . Briefly, Im2shape fits 
each object image with two concentric elliptical Gaussians con- 
volved with a local estimate of the PSF, taking into account 
the background and noise levels. The local PSF is determined 



as the average shape of the three stars closest to each object. 
Shape parameters and ellipticities with their error estimates are 
measured from the deconvolved object images using a Monte 
Carlo Markov Chain (MCMC) optimization technique. 

2.6. Total optical luminosities of the clusters 

There are several ways to estimate the total optical luminosities 
of the clusters from our wide-field images. In all cases, abso- 
lute luminosities include a A;-correction obtained at the clus- 
ter redshift for an ellipti cal-type spectral energy distribution 
dBruzual & Charlotll2003h . 

First, the total luminosity is integrated in circular annuli 
from the "bright galaxy" catalog and corrected for an average 
background luminosity computed in the external areas of the 
images. This is an easy way to compute luminosity and mass- 
to-light ratio profiles, as shown for A 1689 in Paper 1. The main 
difficulty is the determination of the background contamina- 
tion which is estimated from the average density far from the 
cluster The measured values also strongly depend on the con- 
tamination from a few very bright foreground galaxies which 
appear in some cluster fields. To estimate the total luminos- 
ity of the cluster, it is necessary to correct for the magnitude 
cuts of the "bright galaxy" catalog. Assuming that the cluster 
galaxies follow a Schechter luminosity function, the correction 
factor, defined explicitly in Paper 1, ranges from 10 to 30 per- 
cent, depending on the filter and depth of the individual images 
(for example in A 1689 the correction factor is 1.28 in B, 1.11 
in R, 1.27 in /). The "total" luminosity of the clusters is then 
computed in a radius defined by the parameter raoo of the best 
fit model from McAdam (Sect. |4] and Table |5]l. We refer to this 
estimate of the total luminosity as L^' in the rest of the paper. 

We also use the "elliptical galaxy" catalogs, which we ex- 
pect to be cleaner from non-cluster contamination, especially in 
the bright magnitude bins. Another advantage is that these cat- 
alogs are not limited at faint magnitudes and thus include most 
of the elliptical cluster members up to much fainter magnitudes 
than the "bright galaxy" catalogs. Therefore, no incomplete- 
ness factor is necessary in the integration. The correction for 
field contamination, again estimated from the average galaxy 
density outside the cluster and based on the same catalogs, is 
less dominated by bright galaxies than in the previous determi- 
nation. Again, the total luminosity of the ellipticals is computed 
inside the virial radius r2oo and is referred to as L|j" in the rest 
of the paper. 

In practice, the field contamination from a few bright galax- 
ies remains an issue for the total luminosity computed from the 
"bright" galaxies. We therefore use preferentially the luminos- 
ity of the elliptical galaxies of the clusters, which yields cleaner 
light-density maps with a more uniform background distribu- 
tion. We check that the two luminosities do not differ strongly: 
the average ratio L^'/^^" is 1.34, meaning that about 35% of 
the cluster luminosity are missed when only elliptical galaxies 
are considered. Of course, this ratio varies from cluster to clus- 
ter and is also sensitive to differences in the galaxy populations. 
Using L^" only means that we assume that elliptical galaxies 
are best suited for tracing the galaxy content of the clusters and 
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its relation to the dark-matter distribution. Support for this as- 
sumption is provided by several studies showing a strong cor- 
relation between the Ught of early-type galaxies and the mass 
distribution derived from weak-lensing anal yses ( Smail et al] 
1997tlciowe et al.ll998tlGrav et al.l2002l:lGavazzi et alJ20o4 . 



3. Two-dimensional mass maps 

We use an entropy -regularized maximum-likelihood technique, 
implemented in LensEnt2 (Marshall et al. 2002), to obtain 
mass reconstructions for each cluster based on the "faint 
galaxy" catalogs in the R band after PSF correction. Using 
the magnitude-selected "faint galaxy" catalog instead of the 
color-and-magnitude-selected "red galaxy" catalog does not 
change the global shape of the derived mass distribution. To 
avoid building under-dense regions with negative mass den- 
sities a mean background of lOOM0pc"^ is artificially added 
during the reconstruction. Since the mass distribution of clus- 
ters is extended, the individual values of the mass density as 
reconstructed with LensEnt2 are spatially correlated. The ef- 
fective scale of this correlation is controlled by the Intrinsic 
Correlation Function (ICF) of the model which thus sets the 
resolution at which mass structures can be detected. We use an 
ICF of 180" (~ 600/17,5 kpc at the cluster redshift) which rep- 
resents a good compromise between detail and smoothness of 
the mass maps. 

Figs. |4] and |5] show the results of the two-dimensional mass 
reconstructions around each cluster In each plot the signif- 
icance increases by Icr between adjacent contours, with the 
lowest level representing a significance of 2cr above the mean 
background. For each cluster, cr is the average level of the 
noise peaks above the background, obtained from a random- 
ization over the orientation of the "faint galaxy" catalogs used 
in the LensEnt2 mass reconstruction. 200 such random cata- 
logs are produced for each cluster for the statistical analysis 
(see Paper I). The average value of cr for the entire sample is 
CTav - 8OM0pc"2. 

In all cases, the target cluster is detected at high signif- 
icance. Table H] gives the significance of the detection Vpeak 
in units of cr as defined above. Even for the least significant 
detection, A 267, the main peak is detected at about 3.6cr. A 
quantitative assessment of the reality of the mass clumps de- 
tected with LensEnt2 outside the main cluster component and 
at much lower nominal significance would require numerical 
simulations that are beyond the scope of this paper We here use 
the mass maps solely to assess qualitatively the morphology of 
the cluster mass distributions (Table|4]i. We crudely classify the 
clusters as circular or elongated if the mass distribution within 
the ScTav contour has an ellipticity (1 - b/a) smaller or larger 
than 0.2, respectively. 

Also shown in Figs. |4]and|5] are contours of the light density 
of the "elliptical" galaxies as defined in Sect. 12.41 In Table|4]we 
indicate whether there is a good visual correlation between the 
mass and light distributions. Although these are rather qualita- 
tive criteria, most of the "elongated" clusters have both a mass 
and light distribution that is clearly not spherical. 

Globally there is good agreement between the morphologi- 
cal information from our X-ray, strong- and weak-lensing anal- 



yses, as well as the distribution of light in elliptical galaxies. 
We also find good agreement with the overall classification of 
ISmith etalJ toS*): the four clusters A 383, A 963, A 1689 and 
A 1835 correspond to the more relaxed and spherical clusters 
of the sample, while A 2219 and A 68 are also close to this cat- 
egory. 

4. Weak-lensing masses 

4.1. Selection of background galaxies 

One of the main difficulties in obtaining reliable weak-lensing 
mass estimates is to ensure that only background galaxies are 
used in the analysis. The catalogs must be as free as possi- 
ble from contamination from foreground or cluster galaxies to 
avoid attenuation of the weak-lensing signal, averaged in ra- 
dial bins or locally, from galaxies with purely random orien- 
tations. The potential for such attenuation is largest near the 
cluster cores where the galaxy density is highest, causing the 
weak-lensing profiles to be flattened and the total masses to be 
underestimated. For a quantitative test of this effect, we create 
mock catalogs with the same density of sources as observed, 
distributed with random orientations and a Gaussian ellipticity 
distribution similar to the one observed after the PSF correc- 



tion (IBardeau e t al. 2005). These catalogs are then "lensed" by 
a cluster with an NEW mass profile, and average shear profiles 
are built in the same manner as for the observed catalogs (see 
below). The shear profiles are fitted by several mass profiles 
and the total cluster masses estimated in this way compared 
to the input cluster mass. To test the aforementioned attenu- 
ation effect contamination from cluster galaxies is added to 
the lensed catalogs, with a number density profile mimicking 
that of A 1689. The resulting mass profiles and total masses are 
again compared to the initial inputs. We find that even a 10% 
contamination by cluster members in the catalogs can cause 
the reconstructed total mass to be underestimated by up to 20 
to 30%. 

4.2. Fit of the shear profile 

Two parametric mass profiles are used to fit the lensing 
data: a singular isoth ermal sphere (SIS) and an NEW profile 
(INavarro et al. II1997I). The former is characterized by a single 
parameter, the Einstein radius 9e or, equivalently, the velocity 
dispersion cr^,, while the latter is described by two independent 
parameters usually chosen as the scale radius rj and the concen- 
tration parameter c = r2oo/'"s- The virial mass M200 is computed 
as a function of the virial radius r2oo and the critical density pc 
of the Universe at the cluster redshift; 



M20Q = -Trr^oo x200pc(z) 



(1) 



In order to determine the cluster mass distribution from 
the weak-lensing data we fit the s hear pattern using McA dam, 
a Bayesi a n met hod developed by [Marshall et alj ( 120021) and 
Marshal] (l2006f ). McAdam works directly on the PSE-corrected 
faint-galaxy catalogs without any radial binning of the data and 
is consequently more flexible than a fit of a few data points 
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Fig. 4. 15' X 15' CFH12k R-hand images of the eleven clusters of the sample. Thick contours represent the light density of "ellip- 
tical" galaxies selected in the R - I versus R color-magnitude diagram; thin contours represent the mass-density reconstruction 
from LensEnt2 (see text for details). The light contour levels are adjusted for each cluster according to its richness (or central 

Hp-naitv"^ and rantrp. arrmin .^liaht nfFsp.t.'i: hp.twp.p.n thp nn.'iiHnna nf thp Tna.'i:.'i np.aka and thp Inratinna nf thp 






Fig. 5. Continued from Fig. H] 
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Table 4. Global properties of the clusters in the sample: peak surface-mass density in terms of the standard deviation cr of the 
background fluctuations (col. 2); morphological characterization of the mass maps (col. 3); flag assessing the correlation between 
the mass maps and the light maps traced by th e cluster elhpticals (col. 4); X-ray morphology as determined byl Smith et al. (2005) 
(col. 5); "overall classification" assigned by ISmith et al.l (l2005b based on HST lens modeling: for a cluster to be classified as 
"regular" it must contain a single mass component that is centered on the BCG, and the A'-band luminosity of the BCG must 
contribute at least 50% of the cluster's central A'-band luminosity relative to the total cluster luminosity (col. 6). 



Cluster 


Vpeak 


Mass 


M traces L? 


X-ray 


Overall 






morphology 




morphology 


classification 


A 68 


5.5 


Circular 


Y 


Irregular 


Unrelaxed 


A 209 


4.9 


Elongated 


N 


Irregular 


Unrelaxed 


A 267 


3.6 


Elongated 


N 


Elliptical 


Unrelaxed 


A383 


3.6 


Circular 


Y 


Circular 


Relaxed 


A 963 


4.2 


Circular 


Y 


Elliptical 


Relaxed 


A 1689 


9.6 


Circular 


Y 






A 1763 


7.4 


Elongated 


Y 


Irregular 


Unrelaxed 


A 1835 


8.7 


Circular 


Y 


Circular 


Relaxed 


A2218 


6.5 


Elongated 


Y 


Irregular 


Unrelaxed 


A22I9 


5.0 


Circular 


Y 


Irregular 


Unrelaxed 


A 2390 


7.4 


Elongated 


N 







Table 5. Results of the mass-profile fits, obtained with McAdam and using the "red galaxy" catalogs only, for both the NFW and 
the SIS profiles. For all fits the cluster center is assumed to coincide with the location of the brightest cluster galaxy. The Einstein 
radius is computed at the average redshift of the background sources for each cluster (see text for details). M20Q (also referred 
to as the virial mass) is the total mass included in a sphere of radius r2oo. Also listed are the /?-band luminosities of the clusters: 
L^' is integrated within the "bright galaxy" catalog and corrected for incompleteness, while L^" is the total luminosity of the 
"elliptical galaxies" selected within the color-magnitude sequence R- 1 versus R. All luminosities are integrated within the virial 
radius r2oo and error bars are based on the uncertainty of the r2()o determination. 



Cluster 




c 


''200 


M200 




Oe 




^tot 










(hi^ 


Mpc) 


{lO"-h-^MQ) 


(km/s) 


(") 






(10'2 4'Lo) 


A 68 


3.84 


± 1.13 


1.49 


±0.18 


620 ± 


197 


880 ± 65 


12.9 ± 


1.9 


6.6 ± 0.3 


5.5 ± 0.3 


A 209 


3.00 


±0.92 


1.57 


±0.17 


719 ± 


204 


813 ±70 


12.4 ± 


2.1 


8.8 ±0.7 


7.6 ±0.6 


A 267 


4.54 


±2.01 


1.15 


±0.23 


272 ± 


146 


634 ± 116 


6.7 ± 


2.4 


5.3 ±2.0 


4.6 ± 1.7 


A 383 


2.62 


±0.69 


1.32 


±0.17 


419 ± 


146 


619 ±72 


7.4 ± 


1.7 


7.0 ±2.5 


4.9 ± 1.4 


A 963 


8.35 


± 1.25 


1.33 


±0.10 


396 4 


90 


812 ±67 


12.3 ± 


2.0 


5.9 ±0.8 


4.2 ± 0.5 


A 1689 


4.28 


±0.82 


2.25 


±0.14 


1971 4 


336 


1277 ± 37 


31.8 ± 


1.8 


12.8 ±0.5 


11.0 ±0.4 


A 1763 


2.63 


±0.63 


1.93 


±0.14 


1386 4 


263 


932 ± 60 


15.0 ± 


1.9 


12.9 ±0.8 


10.0 ±0.6 


A 1835 


2.58 


±0.48 


2.39 


±0.14 


2707 4 


414 


1240 ±47 


26.6 ± 


2.0 


16.8 ±0.8 


12.1 ±0.9 


A2218 


6.86+1.30 


1.81 


±0.14 


971 ± 


215 


1040 ± 50 


21.2 ± 


2.0 


12.8 ±0.5 


6.8 ±0.3 


A2219 


3.84 


±0.99 


2.25 


±0.18 


2094 4 


435 


1175 ±53 


23.4 ± 


2.1 


14.6 ±0.7 


11.1 ±0.8 


A 2390 


5.26 


± 1.43 


1.74 


±0.17 


943 ± 


246 


1015 ±54 


18.1 ± 


1.9 


9.0± 1.1 


6.6 ±0.8 



obtained by averaging within circular annuli. The output of 
McAdam is a probability distribution of the fitted parameters, 
which is obtained using a maximum-likelihood estimator and 
a MCMC iterative minimization. We use McAdam on the "red 
galaxy" catalogs and fit the mass distribution of each cluster 
with a single component, leaving the mass profile parameters 
(0E for SIS, c and M200 for NFW) as free parameters. A prior 
on the concentration parameter c is included with 2 < c < 10, 
following the results of N-body simulations of cosmological 



structure formation a t the galaxy cluster scale (iBuUock et al 



[2001 ; He nnawi et al.l l2007). As discussed before, any residual 
contamination of the catalog by cluster galaxies, in particular 
in the central area, may reduce the central shear signal and thus 



flatten the deduced mass profile. As no strong-lensing infor- 
mation is included, we add a prior on the lens center by as- 
suming it to coincide with the position of the brightest cluster 
galaxy (BCG). This reduces the number of free parameters in 
the global fitting but is of little consequence for the final value 
of the total mass because the results of the fits are dominated by 
the shear signal at large distance from the cluster center In ad- 
dition to obtaining a global fit, we also compute for each cluster 
the radial shear profile which allows a straightforward assess- 
ment of the strength of the detected shear signal. The eleven 
shear profiles derived from the "red galaxy" catalogs are shown 
in Appendix [A] A detailed cluster by cluster comparison with 
previous mass estimates from independent weak lensing stud- 
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Fig. 6. Mean radial profile of the "red" galaxies, averaged over 
the four "relaxed" clusters (A 383, A 963, A 1689, and A 1835). 
Also plotted is the mean radial profile for the "bright" galaxies, 
representing the cluster members complementing the "faint" 
galaxy sample. As expected the "red" galaxy profile is flat and 
therefore mostly free from cluster contamination. The central 
dip in the faint galaxies profile may be partly due to the magni- 
fication bias acting on the background sources and also to ob- 
scuration effects caused by contaminating bright galaxies (not 
corrected). 

ies is also included in this Appendix. In general, there is good 
agreement between our measurements and previous ones, ex- 
cept in a few cases like A 1689 where there is a wide range of 
measured weak lensing masses, using different methods. The 
most recent ones seem to converge to a value close to the one 
presented in this paper. Note also that seve n out of our eleven 
clusters were imaged bv lDahle et al.l (120021) who fitted their ra- 
dial shear profiles with SIS profiles. Most of the clusters have 
a best-fit value of the velocity dispersion o-jIj much higher than 
the present one, but their analysis was based on shallower data 
taken in a much smaller field of view and so has much larger 
systematic and random errors than our analysis. 

The emphasis of the work presented here is on the total 
mass of the clusters. The resulting best fits as well as the inter- 
nal errors are summarized in Table |5] Given the correlated na- 
ture of the parameters in our analysis and the weakness of some 
of the detections, we focus our discussion on our estimates of 
the virial radius r2oo and the virial mass M2()o which are ro- 
bust results of the weak-lensing analysis. Combining strong- 
and weak-lensing effects increases the accuracy of the mass 
reconstruction close to the center and the fit of the concentra- 



tion parameter c of the NFW profile (iBroadhurst et al. 2005b) 
but requires modifications of the likelihood estimators used by 
McAdam in the central areas. Results from a high-resolution 
study based on this approach will be presented in a forthcoming 
paper (Hudelot et al., in preparation) and will also allow us to 
discuss the distribution of values measured for the concentra- 
tion parameter c across the sample as well as its cosmological 
consequences. 

Fitting the shear profile with the SIS mass profile yields 
a velocity dispersion (Tsbtar which can be converted into an 
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Fig. 7. Projected mass enclosed within a ra dius of 500 / z jp kpc 
derived from the strong-lensing analysis of ISmith et al.l ([2005h 
versus the difference between the weak- and strong-lensing 
masses measured in the same aperture. The weak-lensing 
masses have been converted to an Einstein-de Sitter cosmol- 
ogy with Hq = 50kms ' Mpc"' to allow a direct comparison 
with the strong-lensing masses. 

Einstein radius Oe using the value of the ratio Dls/^^os av- 
eraged over all the sour ces with redshifts est imated from their 
photometric properties (IBardeau et al.ll2005b . 6e does not de- 
pend strongly on the exact value of zs which turns out to be 
Zs - 1-1 for all clusters. We stress that errors due to differ- 
ences in the redshift of the background population from cluster 
to cluster are negligible as we selected the clusters to lie in a 
narrow redshift range so that the scaling factor Dls/^^os varies 
by less than 5% from cluster to cluster. A global error on the 
mean redshift of the sources may introduce a systematic shift 
of the scaling factor and hence the total mass, but its variation 
with (zs) is so small that this cannot account for more than a 
few percent, provided the sources are at redshift at least twice 
the lens redshift. 

A direct comparison between the masses fro m our weak 



lensing analysis and the masses measured by ISmith et al 



(2005) is presented in Fig. [T] The strong-lensing "total mass" 
is in fact the projected mass enclosed within a radius R - 
500 /ig,5 kpc, extrapolated from the strong-lensing modeling, 
whereas the virial mass M200 is measured in a sphere of radius 
r2oo which varies from cluster to cluster In order to account 
for these differences we compute the projected mass inside the 
projected radius R = 500/ijQkpc using the NFW parameters 
determined from the weak-lensing best fit, following the rela- 
tion 



M^^(<R) = 2nr',^, 



Jo ' 



K(y)ydy 



(2) 
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Fig. 8. Comparison between the total mass M200 divided by the 
virial radius r2oo and the velocity dispersion of the fitted SIS 
model. The straight line corresponds to M2oo/'"200 ct^yl 
confirms the correspondence between the two mass estimates 
(see text for details). 



with n = rj mlc. The integral has been c omputed by Bartelmaniil 
(Il996h and lWright & BrainerdI (l2000l) and the projected mass 
can be rewritten in terms of the virial mass M200 and the con- 
centration parameter c as 



M200 IR 
ln(l+c)-c/(l+c) \h 



(3) 



with Tix) defined as 



vr 



: artanh 



\ + X 



y Vx2- 1 



: arctan 



X- 1 
x+\ 



(i) 



(^< 1) 

{x^D 

{x>\) 



(4) 



The comparison between the two masses is shown in Fig.|7] 
and shows good agreement. No obvious bias appears in the 
sample although there seems to be a dichotomy between two 
families of clusters. A 209, A 1763 and A 2219 have weak- 
lensing masses that appear to be overestimated in comparison 
with those of the other clusters. By contrast, the three clusters 
identified as "relaxed" clusters both in the weak- and in the 
strong-lensing analysis (A 383, A 963 and A 1835) show good 
agreement between the two mass measurements. 

In order to check the validity of our mass estimates, we 
also compared masses deduced from the radial fits of the two 
different parametric models (NFW and SIS). The total mass 



enclosed within the virial radius is M200 for the NFW profile 
and is proportional to o\^^ X ^200 for the SIS profile. Fig. [8] 
compares these two quantities which show on average the same 
behavior, so we are confident that using either the virial mass or 
the mass deduced from the isothermal profile does not change 
the validity of the scaling relations discussed below. Note that 
in both cases we use the 3D total mass enclosed within the 
virial radius r20(). 

5. Global correlations 

With all measurements in hand we now investigate the correla- 
tions between the lensing mass and immediate cluster observ- 
ables like the optical luminosity or X-ray characteristics such 
as luminosity and temperature. 

5.1. Correlations between mass and optical cluster 
properties 

Mass and luminosity are strongly correlated quantities as can 
be seen in Fig.|9] The M/L ratio is usually representative of the 
dynamical state of a cluster and is a tracer of its star-formation 
history. One of the current issues of cluster research is whether 
the M/L ratio is a constant and universal value at least for rich 
clusters, or w hether it increases w ith mass, as recently sug- 
gested by, e.g.. |Popesso et al.l (l2005b . 

We here compute the M/L ratio as the virial mass within 
the virial radius r2oo (^200) divided by the luminosity L^" as 
defined above. We use the virial mass instead of the 2D pro- 
jected mass MjJJ) because the latter depends on the concentra- 
tion parameter c, which is poorly constrained in our study. We, 
however, convince ourselves that changing c from 2 to 10 de- 
creases the deduced projected mass by only less than 20%. The 
optical cluster luminosities are carefully corrected for back- 
ground contamination as well as Galactic extinction and are 
also ^-corrected. We do not apply any correction for luminos- 
ity evolution which is likely to be small at the low redshift of 
our target clusters. To obtain a global M/L ratio for all clusters 
we apply the average correction factor to the total luminosity, 
1 .34. In practice, we find an average value for the whole cluster 
sample before correction for the total luminosity 



<M/4') = (133±50)/i7()(M/L)c 
or globally for all galaxies 
<M/4°') = (100±38)/!7o(M/L)c 



(5) 



(6) 



These values are close to the value o f M/Ly - 



180+^[0 h (M/L)q obtained by lSmail etalJ (11997!) in their early 



weak- lensing study o f a sam ple of rich clusters. On the other 
hand, ICarlberg et al. ( 19961 119971) find an average value of 



(M/Lr) = (289 + 50) hm (M/L),© for the M/L ratio of rich 
clusters from the CNOC survey (a sample of 14 clusters span- 
ning a redshift range [0. 17 - 0.54]). Theirs are, however, global 
dynamical values which overestimate the M/L ratio measured 
at the virial radius r2oo by approximately 20% ( Carlberg et al] 
19971) . Although, all considered, their average M/L ratio is still 
higher than the one found in our analysis, the discrepancy is 
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Lr-^" (10'2 Lo) 

Fig. 9. Mass versus optical luminosity for the clusters in our 
sample. The mass is computed at the virial radius raoo derived 
from the best weak-lensing fits obtained with McAdam. The op- 
tical luminosity is computed in the R band for the "elliptical" 
galaxies only, selected from the color-magnitude relation, and 
has been corrected for background contamination. Both quan- 
tities are projected (i.e. integrated along the line of sight) and 
are therefore comparable. The dashed line represents a constant 
M/L ratio of 133 in solar units while the solid line represents 
the best-fit power law M cc L' 

barely significant when all sources of error are taken into ac- 
count. For A 2390, the only galaxy cluster in common between 
the two samples, the velocity dispersion as well as the M/L 
point toward a higher mass than the weak lensing mass mea- 
sured in this paper However, their determination of the virial 
radius differs significantly from ours. As a consequence, the 
virial mass and the total luminosity are not measured in the 
same area and the comparison is not conclusive. 

If we drop the assumption of a constant M/L ratio for clus- 
ters, Fig.|9]shows a strong correlation between mass and light, 
steeper than for a constant M/L, with massive clusters having 
higher M/L ratios. A log-log fit of the mass-luminosity rela- 
tion, including errors on the mass measurements shows a power 
law dependence with an index a = 1.80 + 0.24, or equivalently 



M/L oc l' 



0.80±0.24 



(7) 



The tendency of increasing M/L ratio with virial mass is sig- 
nificant, even though the mass range of the present sample 
is rather limited. Previous analyses have already shown simi- 
lar trends, in particular in a detailed comparison between dy- 
namical mass estimates from the SDSS and luminosities com- 



ing fr om either optical or X-ray measurements (IPopesso et al 



2005h . Their sample spans a wide mass range, from groups 
up to massive clusters, and although there is a large disper- 
sion of their data points, the correlation between M/L and 



mass is confirmed, but with a lower slope than in the present 
study. However, although the mass range of their sample is 
wider than ours their dynamical mass determination may not 
be as accurate as the weak lensing masses in the present anal- 
ysis. Complementary lensing mass measurements for lo wer 
mass clusters have been obtained by ^ Parker et al. (2005) us- 
ing the CNOC galaxy groups. They show some evidence for 
an increase in the M/L ratio from poor to rich galaxy groups. 
This result is in agreement with theoretical predictions ob- 
tained by comparing the Press-Schechter mass function with 
the observed luminosity function dMarinoni & Hudsonll20()2l) . 
The M/L ratio in the high mass range, typical of rich clus- 
ters of galaxies, scales as iP-^'^^-^^^ in close agreeme nt with our 
observations. However. [Marinoni & Hudso assumed a 

rather high value of erg - 0.9. Following the latest WMAP3 
analysis, a lower erg = 0.75 might change their M-Lo^t relation 
by lowering the normaUzation but not the power-law exponent 
(lReiprichll2006l) . 

This scaling behavior is predicted in semi-analytical mod- 
els of galaxy formation and has been interpreted as a de- 
crease in galaxy formation efficiency in rich a nd dense environ 



ments due to the long cooling time of hot gas dKauifmann et al 



19991). Indeed, recent X-ray observations have shown that in 
the centers of dense clusters the gas cools less efficie ntly than 



predi c ted by the standard cooling flow model ( Peterson et al 



20011 iFabian et al.ll2003b . By injecting energy in the form of 
radio jets, AGN may be responsible for switching off" cooling 
at the centers of massive haloes, thus preven ting the formation 
of very bright structures ( Bower et al.ll2006l) . Also, star forma- 
tion is quenched in newly infalling galaxies through a variety 
of physical processes, some of which, e.g. ram pressure strip- 
ping of_gaSj_arefliemore efficient the denser the environment 
(e.g. lTreuetaTll2003h . 



5.2. Correlation between optical and X-ray luminosities 

There is some correlation between the clusters' optical and X- 
ray luminosities, although the luminosity range is not very ex- 



tended in our sample of bright X-ray clusters. iPopesso et al 



(l2005l) performed linear fits between these quantities using a 
large sample of low-redshift groups and clusters of galaxies 
from the RASS-SDSS Survey over more than two decades in 
luminosities. For the correlation between optical and X-ray lu- 
minosities, they find a best fit law; 

logCLopt/lO'^Lo) = 0.64 log(Lx/10'*'^ ergs"') + 0.45 + 0.15 . 

Instead of doing similar fits with our sample, which may have 
little physical meaning, we simply compare the averaged opti- 
cal luminosity (L'°' = 10.2 x IO'^Lq) at the mean X-ray lumi- 
nosity of the sample ((Lx) = 10.5 x lO"''* erg s"'). This value 
corresponds c losely to the value predicted by the relation of 
Popesso et alJ (E005). We therefore confirm their normalization 
from an independent sample with great accuracy and demon- 
strate that this normalization does not change out to z ~ 0.2 at 
high Lx- 
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Fig. 10. Left: Weak lensing 3D virial mass M200 versus X-ray luminosity. The best fit line has a slope a - 1.20 ± 0.16 and is 
discussed in the text. Right: Weak lensing 3D virial mass M200 versus X-ray temperature. The straight line corresponds to a 
M200 °^ T^^^ relation while the dashed line corresponds to the best fit power law relation M oc 7'4 6±o.7 T he virial mass M?n n 
is derived from the best weak len sing fits obtain ed with McAdam, temperatures are derived from XMM data (IZhang et al.lbOOTF) . 
completed by ASCA data from dOta & Mitsudai 120041) for A 2219. In both plots, the 4 clusters with cooling core or relaxed 
properties are marked wit empty boxes. 



5.3. Correlation between cluster mass and X-ray 
luminosity 

We now compare the deprojected mass derived from the weak 
lensing analysis with the cluster X-ray properties. We consider 
the virial mass M200 enclosed within rioo, the result of the 
weak lensing best fit performed with McAdam. Bolometric X- 
ray luminosities are measured from XMM-Newton data and 
integrated up to 2.5r5oo, a radius w hich is close to twice our 
virial radius ^200 dZhang et al .l2007h . They are listed in Table[T] 
Fig. [To] shows the correlations between M200 and the cluster 
X-ray properties. A linear regression between mass and lumi- 
nosity (including errors in both axes in log-log space) gives the 
scaling relation 



1.20±0.16 



(8) 



M20(,/10'2Mo =40!?3 (Lx,boi/1044ergs-') 

As the slope is close to 1, one may also consider a constant 
mass to (X-ray) light ratio for which we obtain an average value 
of M/Lx = 67±24in units of (lO'^Mo/lO'^'^erg s"'). Inverting 
the Lx - M relation gives 



ix,boi/10''*ergs-' 



2.16!°:^^(M2„o/10'4Mo)' 



0.83±0.11 



(9) 



which seems in contradiction with some previous determina- 
tions of the correlation between mass and X-ray luminosity 
where the slope of the Lx - M relation is on the order of 



1.3 to 2.0 dReiprich&Bohringeri 120021: iPopesso etalj 12005 



MaughanI 120070 . From the theoretical point of view, a rela- 
tion Lboi M'^^^ is expected in a self-similar model for clus- 
ters, while including "pre-heating" in the physical processes 



of cluster formation lea ds to a relation given by Lboi M"^^ 
(lEvrard & Henrylll99lh . An accurate determination of the ob- 
served Lx - M relation is quite difficult with our present mass 
estimates, also because the perturbation induced on the mea- 
surements by the central cooling core of some clusters modifies 
the sample properties. Note also that our present normalization 
of the L - M relation is only 60% higher than the one deter- 
mi ned by Maughan ( 2007), contrary to the normalization found 
bv^Popess o et al.i (i2005i) which is lower than our value by a fac- 
tor of 3 at IO'^Mq but their L-M relation presents a very differ- 
ent slope, close to 2 when using bolometric X-ray luminosities. 
Going further in this analysis is presently difficult because of 
the limited size and the small mass range of our sample. In ad- 
dition, our sample is X-ray flux-limited which may introduce 
some bias into the relations involving Lx, because for a given 
mass, only the high luminosity clusters are selected (.ReiprichI 
20061). 



5.4. The mass-temperature relation 



The X-ray tem peratures have be en measured on detailed XMM 
observations tZhang et al ] |2007h . The global values are the re- 
sults of the volume averaged radial temperature profiles be- 
tween 0.2 - 0.5r2oo- As shown in Fig. [10] we presently have 
too large uncertainties and dispersion in the measurements to 
correctly fit a power law for the M - T relation. Instead, we 



S. Bardeau et al.: Weak lensing survey of X-ray luminous galaxy clusters. II. 



15 



simply fix the slope to what is expected if clusters are in hydro- 
static equilibrium, M oc T^l-. Then we find a normalization 



M20(,/10'^Mo = 0.44!|;:^^ 



llkeV/ 



(10) 



This normalization deserves some comment because it is one 
of the few attempts at fixing i t with weak lensing mass mea- 
surements. Pedersen & Dahle ( 20061) also provided a tentative 
determination of the M -T relation with weak lensing masses. 
From their whole sample of rich clusters, spread over a large 
redshift range, they find 

Mskev = (0.78 + 0.14) X 10'^ /i"' Mq . 

The normalization translates to a value of 0.49 ± 0.09 in our 
units (Eq. [TOl l and is in very good agreement with the present 
work. Our value is also quite close to a recent prediction ob- 
tained from a sample of simulated clusters realized with the 
most recent imp lementation of the Tree-SPH c ode GADGET2 
('SDringef2005'). For the most massive clusters, Ascasibar et al 
(,2006.) find a M - r relation of 



Mzoo/lO'^^Mo = 0.473 /17 



3/2 



,lkeV/ 



(11) 



also in good agreement with our result, although their relation 
is established at z = 0. However, the authors claim that their 
normalization is lower than the values found in previous sim- 
ulations and explain this partly by improvements of the treat- 
ment of entropy conservation. 

5.5. Discussion 

The scatter in the M-T and L-T relations should be represen- 
tative of the diversity in the cluster histories but it is not easy 
to analyze because our sample is small. The individual mass 
measurements still have large uncertainties associated to the 
weak lensing method: masses are underestimated if the cata- 
logs are not completely cleaned from foreground or cluster con- 
tamination and the low density of background sources used for 
the weak lensing reconstruction adds another source of noise 
which strongly limits the accuracy of the measurements. These 
systematic uncertainties, which are difficult to quantify without 
detailed simulations of mock catalogs, are presently not taken 
into account in this analysis . The remaining uncertainty on the 
concentration parameter adds at least another 20% uncertainty 
on the total mass, included in our present mass errors budget. 
All in all, weak lensing masses cannot presently be determined 
with accuracy better than 40 to 50%. However, the tight corre- 
lation between mass and optical luminosity suggests that some 
of these biases partly cancel when measuring the mass or the 
total luminosity; this is the case for all geometric departures 
from spherical symmetry (ellipticity of the light/mass distri- 
butions or projection effects). Systematic biases in the mass 
determinations should not change the slope of the scaling re- 
lations dramatically while some uncertainty in their normaliza- 
tion remains. But since the slopes of the Lx - M, Tx - M and 
Lop, - M relations are all shallower compared to previous re- 
sults, we may also suspect some scale dependant biases, most 



probably at the low mass end of our sample, where the weak 
lensing measures are the most difficult to caracterize. Further 
improvements in the weak lensing methodology are in progress 
and may help clarifying this possible bias. 

Moreover, although the sample was initially selected for 
its homogeneity, at least in X-ray properties (Lx > 4 x 
10"^"* erg s"'), at least three out of the eleven clusters 
(i.e. about 25% of the sample) present a strong central cooling 
core (namely A 383, A 963 and A 1835) which perturbs the to- 
tal X-ray luminosity (clusters are overluminous) and the X-ray 
temperature (clusters are too cool) for a g iven mass. This ef- 
fect is partly taken into account in the wayl Zhang et al. (2007) 
measured the X-ray temperatures, excluding the X-ray signal 
in the central core (r < 0.2rgoo). Another five clusters show 
signs of non- sphericity in the mass distribution. Following 
Popesso et alj ([2005) we suspect that most of the dispersion in 
the scaling relations is due to the intrinsic dispersion on the X- 
ray properties of these rich clusters although the present mass 
uncertainties add a significant fraction of the total scatter. For 
the optical luminosity of clusters the link with the total mass is 
tighter This is an indication that early-type galaxies are good 
tracers of the mass in clusters and that the dispersion in M with 
Tx is real. 

6. Conclusions 

We have presented the first weak-lensing analysis of a homo- 
geneously selected sample of eleven X-ray luminous clusters 
in a narrow redshift slice at z ~ 0.2. Using wide-field imag- 
ing in three bands (B, R and /) covering up to 5 Mpc in radius 
around the cluster targets we are able to disentangle between 
foreground, cluster and background galaxies. The weak lens- 
ing signal is always well detected up to 2 Mpc, but generally 
extends out to the edge of the field at lower signal-to-noise. 

The weak lensing methodology used in the paper has al- 
ready been tested and validated on the cluster A 1689 which 
yields the highest s ignal-to-noise detection in our sample 
(IBardeau et al.l 120051) . To measure galaxy shapes and correct 
them for PSF anisotropy and circularization, we used the 
1m2shape tool. Our reliance on 1m2shape is j ustified by the 



results of the Shear Testing Program (STEP, Hevmans et al 



2006h which fi nds it to be a promi sing alternative to the popular 
KSB method (Kaiser et al. 1995). Cluster masses were com- 



puted using the McAdam software (lMarshallll2006 ) which per- 
forms a two-dimensional fit of the individual galaxy shape in- 
formation of the "background" galaxy catalog. We were thus 
able to determine: 1) M200 derived from fitting NFW profiles 
(note that only weak constraints could be placed on the concen- 
tration parameter, which is generally degenerate without any 
additional strong lensing constraints); 2) the velocity dispersion 
cry derived from fitting SIS profiles. Although the details of 
the mass profiles are not well constrained by the weak-lensing 
analysis, in particular close to the center, the total mass is ro- 
bust with respect to variations in the concentration parameter 
of the NFW profile. 

We believe that a number of limitations are still affecting 
our lensing mass estimates: 1) we have only limited knowledge 
of the redshift distribution of the galaxies used to measured the 
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shear signal - this may plague the measurement of the shear 
signal towards the cluster center, although the narrow redshift 
range of the sample (Az/z ~ 10%) and the homogeneity of the 
weak lensing data make a direct comparison between the dif- 
ferent clusters possible; 2) we could recover the shear signal at 
best to the 20% level due to a relatively small usable number 
density of background galaxies (about 10 galaxies per square 
arc-minute after taking into account a color selection) and a 
probably not perfect PSF correction; 3) our lensing mass deter- 
mination assumes circular symmetry which may in some cases 
be a poor representation of the cluster morphology. Fortunately, 
thanks to the large spatial extent of the weak-lensing detection 
the mass-sheet degeneracy in our mass measurements is mini- 
mized. 

As for the total masses M200 of the clusters in our sample, 
we investigated the relations between mass. X-ray and optical 
observables. Unlike most previous attempts at calibrating these 
relations, we use masses measured directly from their gravita- 
tional effects employing methods and data that are completely 
independent of the X-ray measurements. In particular, the er- 
rors in the X-ray observables and in the weak-lensing masses 
are uncorrected. We summarize the main conclusions and re- 
sults of the study: 



- The optical M/L ratio presents a strong correlation with the 
cluster luminosity, with MjL oc /^o so±o.24 most mas- 
sive and luminous clusters thus have the highest M/L ratio. 
This reflects a change in galaxy formation efficiency in rich 
clusters. 

- There is a strong correlation between mass and X-ray lumi- 
nosity with Lx MjQQ^^^'' ". The logarithmic slope is sig- 
nificantly smaller than found in previous attempts to com- 
pare both quantities. A better understanding of the behavior 
of this relation is crucial in view of future large X-ray sur- 
veys of clusters for most of which only X-ray luminosities 
will be available. 

- The mass range of our cluster sample is too small to cor- 
rectly fit a M - r relation so we fix the logarithmic slope 
to 3/2 and concentrate on determining the normalization 



to find Mzoo/lO'^Mo = 0.473 /i^^ (Fx/Ike V)^^l This nor- 
malization is very close to the value predicted from numeri- 
cal cosmological simulations of cluster formation and evo- 
lution as well as the observed normalization from X-ray 
measurements ( Ascasibar et al ] |2006HSpringell2005h . This 
good agreement also suggests that evolutionary effects are 
negligible between z - and z - 0.2. 
- The scatter in the M - T relation is still large and difficult 
to disentangle from uncertainties in the measured masses 
introduced by limitations in the current lensing analysis. 
Analysis of larger samples of clusters with, ideally, better 
weak-lensing data and comparison with simulated data will 
be required to conclusively address the impact of the hydro- 
dynamical state of clusters and the reliability of X-ray mea- 
surements. It is, however, akeady clear from this study that 
our eleven clusters differ significantly in terms of global 
morphology, dynamical state, mass concentration, and thus 
possibly merging histories. 



In the near future, progress on this kind of analysis may 
be achieved in different ways. Better constraints on mass pro- 
files and particularly on the concentration c will be discussed 
in a forthcoming combined analysis of strong- and weak- 
lensing data where the number and positions of strongly lensed 
or multiply imaged galaxies unambiguously determines the 
slope of the mass profile in the cluster cores (Hudelot et al., 
in preparati on). First results hav e already bee n obtained in 



Abefl 1689 (iLimousin et al.ll2007l) and Abell68 (Richard et al 



2007). Analysis of larger samples of clusters with current facil- 
ities may provide clues on cluster physics by statistically re- 
ducing the scatter in the various scaling relations. However, 
the final limitation may arise from the intrinsic scatter of the 
X-ray properties (non-equilibrium processes, cluster mergers 
which tend to increase both X-ray temperature and luminosity) 
or the lack of accuracy in our lensing measurements. In con- 
clusion, better lensing measurements will likely be obtained, in 
the near future, with very deep multi-color imaging using wide 
field cameras (Subaru/SuprimeCam or CFHT/Megacam) or in 
the more distant future using a space-based wide field imager 
such as the SNAP telescope. 
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Appendix A: Shear profiles of the cluster sample 



The shear profiles for all clusters are shown in Fig. lA.ll They 
are based on the "red galaxy" catalog and feature statistically 
independent data points. We also plot the results of the McAdam 
fits obtained with the two mass profiles discussed in the main 
body of this paper (SIS and NEW, see Table |5]). The shear sig- 
nal in the central bin is generally similar to or weaker than the 
signal in the second bin. This is likely due to residual contami- 
nation from cluster members or foreground galaxies in the cat- 
alogs. 



central part of the mass model. The strong-lensing model and 
the weak-lensing signal appear to be in good agreement in this 
cluster although they are in disagreement with the previously 
measured weak lensing mass of ,Dahle et al.. (,200Z) . 



B.2. Abell209 

Abell209 is dominated by a bright cD galaxy that is elon- 
gated in the NW-S E direction. There are no obvious giant arc 
systems, although iDahle et alJ (l2002h mention an arc candi- 
date embedded in the envelope of the central galaxy. The in- 
ternal dynamic s of th e cluster has been studied in detail by 
Mercurio et al.l (l2003b who find a high value of the line-of- 
sight velocity dispersion (cr\os ^ 1400 km s"'). The presence 
of a velocity gradient along the main extension of the dis- 
tribution of the cluster galaxies and evidence for substruc- 
ture and dynamical segregation suggest that we are observing 
this system in a late merger phase. This inte rpretation is sup 



porte d by the recent detection of a radio halo (iGiovannini et al 



20061) which implies that this cluster is a dynamically immature 
merger An indepe ndent weak-lensing analysis of Abell209 
was performed by iPaulin-Henriksson et al.l (l2007h using the 
same data retrieved from the CFHT archives. They measured 
a virial mass M2(K) = 7.7+4-^ X 10^* Mq inside a vmal radius 
'"200 = (1.8 ± 0.3) Mpc, values that are in good agreement with 
our own measurements although they were obtained with a dif- 
ferent weak-lensing pipeline and me thodology. They a re also 
compatible with previous estimates dDahle et alj|2002h . Their 
results also show a strong elongation of the cluster. 



B.3. Abell267 

Abell267 is very similar in appearance to Abell 209. The clus- 
ter is dominated by a giant cD galaxy with large ellipticity. 
There are no giant arcs and no obviously lensed background 
galaxies (Smith et al. 2005). The weak-lensing signal is the 
weakest of the sample, making this one of the least massive 
systems studied here. 



Appendix B: Individual properties of the clusters 

B.1. Abell 68 

Abell 68 contains a large cD galaxy elongated in the NW-SE 
direction. About 1 ' to the North- West of the cD is a compact 
group of about five bright galaxies, which gives the system a 
bimodal appearan ce. Many blu e arcle ts can be seen around the 
cluster center and I Smith et alJ (120051) provide several redshift 
identifications for these objects. Most notable is a triple arc lo- 
cated just east of the center of the cD galaxy. Recent integral- 
field spectroscopic observations revealed very strong Lyman- 
a emission correspo nding to an extended object at z - 2.63 
dRichard et al.ll2007b . This dedicated spectroscopic survey of 
all the lensed galaxies in this cluster also provides a catalog 
of 27 images with spectroscopic identification and redshifts 
ranging from 0.4 up to 5.4, which considerably improves the 



B.4. Abell 383 

Abell 383 is dominated by a nearly circular cD galaxy and 
shows a rich and complex system o f giant arcs and arclets 
dSmith et"aL 2001 : Smith et al. 2005 ). These arcs were used 
to reconstruct the central mass distribution of the cluster with 
high accuracy. With five multiple-image systems identified in 
the cluster center and three more faint, tentative systems, the 
mass distribution in the core is tightly constrained. The central 
slope of the mass profile is surprisingly steep, but is consistent 
with the excess mass due to the strong cooling core which feeds 
baryonic mass to the cD galaxy. The weak-lensing signal, on 
the other hand, does not present the same characteristics, lead- 
ing us to conclude that this cluster is highly concentrated but 
not widely extended. 
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Fig. A.l. Shear profiles of the 1 1 clusters, using the "red galaxy" catalog. The results of the fits with two different mass profiles 
(SIS and NFW) using McAdam are shown. 



B.5. Abel 1 963 

The center of Abell 963 is dominated by a cD galaxy and con- 
tains two gi ant arcs to the North a n d the South of th e cD, 



Ellis et al 



(119911) mea- 
0.771, while the 



respectively jLaverv & Henrvlll988l) 
sured a redshift for the northern arc of z 
southern arc, of very blue color, has not been ide ntified spec 
trosco pically. From their strong-lensing model, ISmith et al 



20051) argue that the southern arc is a group o f singly imaged 
galaxi es rather than a multiple-image system. Laverv & Henry 



(Il988h measured a velocity dispersion of ISSO^I^j^Skm s ^ from 



36 cluster members. The weak-lensing signal is remarkably 
regular and consistent with an NFW profile even at small radii. 
The spatial distribution of the cluster members is highly circu- 
lar and characteristic of a well relaxed cluster 

B.6. Abell1689 

Abell 1689 is a well known gravitational lens. It is a very 
rich and luminous cluster, dominated by a compact group of 
bright galaxies and a central giant elliptical (gE) galaxy. Deep 
HST/ACS images reveal a strong over-density of faint com- 
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pact galaxies in the periphery of this galaxy (IMieske et al 
20041) . highlighting its dominance at the center of the clus- 
ter. A dynam ical study of this cluster was performed by 
Teague et al.l ( 1990) who found an extremely high velocity dis- 
persion (cTios = 198 9km s ~ ') pro bably arising from a com- 
plex merger iGirardi et al.l (1997) confirmed this initial as- 
sumption with a detailed study of substructure in the clus- 
ter and showed that the measured velocity dispersion proba- 
bly contains a systematic peculiar velocity component caused 
by an ongoing cluster merger. X-ray properties, such as the 
low gas mass fraction in this cluster , confirm this inter- 
pretation (lAndersson & Madejskil l2004l) . Gravitational lens- 
ing features have b een st udied extensively in Abell 1689. 



Clowe & Schneided ( l200lb and iKing et alj (l2002h present 



weak-shear measurements in the cluster which give a global 
mass profile close to an NFW profile, compatible with the X- 
ray data, but with a total mass lower than our present measure- 
ment by nearly a factor 2. The most spectac ular results come 



from the deep HST/ACS images obtained by iBroadhurst et al 



(l2005al) which display a number of arcs and arclets, making 
this cluster the most spectacular cluster lens with one of the 
largest Einstein radii observed in clusters. iHalkola et al. I (l2006h 
also built a non-parametric strong-lensing mass model and in- 
cluded external weak-lensing constraints to derive a total virial 
mass, with a value that is 50% higher than our measurement. 
More recently, an extensive spectroscopic survey of the arcs 
in this cluster has been conducted (Richard et al. 2007b, in 
preparation) yielding spectroscopic redshifts for about 2/3 of 
the 32 multi ple image syste r ns ide ntified in the cluster Using 



these data, Limousin et al 



(120071) constructed an improved 
mass model combining constraints from both strong and weak 
lensing. They found a value of 7.6 +1.6 for the concentra- 
tion paramet er c, similar to what is expected from numerical 
simulations ( Bullock et al.l EoOl). a total virial mass M200 = 
(1.32+O.2)xlO'^M0, and a virial radius r2oo = (2.16±0.1)M/:>c, 
in excellent agreement with our results. 



B.7. Abell1763 

Abell 1763 has a central cD galaxy but otherwise the clus- 
ter center is comparatively ill defined, with "chains" of bright 
galaxies heading off in at least three directions. There are no 
obvious gravitational-arc systems in the cluster center The 
weak-lensing signal is rather weak even close to the center re- 
sulting in a low value for the total mass. The two-dimensional 
mass reconstruction clearly shows signs of bimodality in this 
cluster, with the main component centered on the cD galaxy 
and a second component 4' to the west . This second mass struc- 
ture was not detected by IPahle et al. (2002) because it is be- 
yond their maximum radius of investigation. 

B.8. Abell1835 

Abell 1835 is dominated by a giant elliptical galaxy slightly 
elongated in the North-South direction. The galaxy distribution 
is regular and the global appearance is that of a well relaxed 
cluster. A 1835 is the most X-ray luminous cluster in the BCS 



catalog dEbeling et al. I I2OOOI) and therefore of the present sam- 
ple. From the Lx - Tx relation it should also be the hottest 
cluster; however, its measured X-ray temperature is relatively 
modest ( Table ^ and dr ops to about 4keV in the cooling-core 
region dMcNamara et al.ll2006l) . From a lensing point of view, 
A 1835 is a well known strong-gravitational lens, with many 
very thin long gravitational arcs seen in the HST/WFPC2 im- 
age. Several of them have been iden t ified spectroscopic ally 



dRichard et alj 120031: IPello et alj |2004 ISmith et all I2OO6I) . m 
particular a controversial z = 10.0 gala xy. A weak lensing 
study of this cluster was carried out by Clowe & Schneiden 
(l2002h with the Wide Field Imager (WFI) on the ESO/MPG 
2.2 m telescope. The isothermal fit agrees well with our work, 
as does the NFW r2oo value which gives a total mass of the 
cluster in very close agreement with our measurement. 



B.9. Abell 221 8 

Abell 22 18 is arguably one of the most famous cluster lenses 
with an extraordinary number of arcs and arclet s in its cen- 



ter. A lens model for this cluster was presented by lKneib et al 



(1 19961) and required a bimodal central mass distribution, with 
one mass component centered on the cD galaxy and a sec- 
ond one center ed on a brigh t galax y about 1.5' to the south- 
east of the cD. iGirardi et al.l ( 119971) analyzed the distribution 
of 50 galaxy redshifts in A 22 18 and found evidence for two 
groups of galaxies superimposed along the line of sight, which 



they i dentify with the mass clumps modeled by iKneib et al 



(Il996l) . Deep HST/ACS images revealed the nature of several 
strongly lensed galaxies and in particular a z ~ 7 g alaxy candi 



date, r ecently confirmed by Spitzer observations (IKneib et al 



2004t lEgami et al.l 120051) . The weak-lensing signal detected 



in this cluster is one of the strongest in our sample, and the 
two-dimensional mass distribution is perfectly matched to the 
galaxy distribution, with an ellipticity and orientation which 
reflects the central bimodality of the mass distribution. 



B.10. Abell 221 9 

The optical morphology of Abell 2219 is remarkably simi- 
lar to that of Abell 2218, with a dominant cD galaxy and a 
second bright elliptical galaxy at ~ 1' to the South- West. 
A number of gravitational arcs can be seen, most notably a 
straight arc between the two brightest cluster galaxies and a 
very thin and elongated arc to the North- West. The presence 
of two mass clumps as well as the elongated X-ray distri- 
bution indicate a non-relaxed cluster in the process of merg- 
ing of several substructures. A gravitational depletion s i gnal a t 
near-infrared wavelengths was detected by Gray et al. I i200flh . 
The mass deduced from this effect was fitted with a singular 
isothermal mass distribution with cr ~ 800 km s"' and is con- 
sistent with our measurement. However, our measurement is 



slightly lower than the one obtained by lDahle et alJ (l2002h . The 
two-dimensional mass reconstruction reflects the general elon- 
gation of the cluster, although this is less significant than in 
Abell 22 18, for example. 
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B.11. Abell2390 



Abell2390 is a cD-dominated cluster with several arcs in its 
center A chain of fairly bright galaxies extends to the North- 
West of the cD galaxy. The arcs on this side of the cD are 
straight (Pello et al. 1991) and confirm the extension of the un- 
derlying mass distribution in this direction. About 3' to the 
East of the cluster center lies an extended group of galax- 
ies. This cluster has a surprisingly high velocity dispersio n 
dLe Borgne et all Il99ll Fee et al.lll99a iBorgani etaLfll999l) . 
X-ray observations reve al a strongly elongated gas distribution 
along the NW-SE axis dPierre et al.| [l996; Allen et al. 2001), 
consistent with the optical galaxy distribution. Our own two- 
dimensional mass reconstruction has a similar elongation in 
the same direction supporting the notion that A 2390 is a non- 
relaxed cluster, accreting one or two mass clumps along the 
same axis. 



